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ABSTRACT 

We present the study of a large sample of early-type dwarf galaxies in the Coma 
cluster observed with DEIMOS on the Keck II to determine their internal velocity 
dispersion. We focus on a subsample of 41 member dwarf elliptical galaxies for 
which the velocity dispersion can be reliably measured, 26 of which were studied 
for the first time. The magnitude range of our sample is —21 < Mji < — 15 mag. 

This paper (paper I) focuses on the measurement of the velocity dispersion 
and their error estimates. The measurements were performed using pPXF (pe- 
nalised PiXel Fitting) and using the Calcium triplet absorption lines. We use 
Monte Carlo bootstrapping to study various sources of uncertainty in our mea- 
surements, namely statistical uncertainty, template mismatch and other systemat- 
ics. We find that the main source of uncertainty is the template mismatch effect 
which is reduced by using templates with a range of spectral types. 

Combining our measurements with those from the literature, we study the 
Faber-Jackson relation (L oc cr") and find that the slope of the relation is 
a = 1.99 ± 0.14 for galaxies brighter than Mn 2± — 16 mag. A comprehensive 
analysis of the results combined with the photometric properties of these galaxies 
is reported in paper II. 

Key words: galaxies: clusters: individual: Coma; galaxies: elliptical and lenticu- 
lar, cD; galaxies: dwarf; galaxies: kinematics and dynamics; galaxies: fundamen- 
tal parameters; galaxies: evolution 



* Based in part on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Institute, 
which is operated by the Association of Universities for Research in 
Astronomy, Inc. under NASA contract NAS 5-26555. These obser- 
vations are associated with programme GO 10861. 
t Some of the data presented herein were obtained at the W.M. 
Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of Cal- 
ifornia and the National Aeronautics and Space Administration. The 



1 INTRODUCTION 

A clear understanding of dwarf galaxies and their relation 
to their most massive counterparts is essential for testing 
galaxy formation models. Although the lowest luminosity 
dwarf galaxies (AI^ ~ — 9) appear to contain large amounts 
of dark matter (Aaronson et al. 1983; Mateo 1998) for 
brighter dwarfs the situation is less clear. De Rijcke et al. 
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(2006) find that the dwarf companions to M31, at M„ ~ 
—14, contain around 40-50% dark matter by mass within 
the inner two half-light radii. Geha et al (2002) find that six 
Virgo cluster dwarfs at Mv ~ —16 do not show evidence for 
dark matter within an effective radius. Toloba et al. (2011) 
observe a larger sample of Virgo dwarfs, and again find that 
they are not dark matter dominated. Here we observe a sam- 
ple of Coma cluster dwarfs in a somewhat brighter magni- 
tude range (21 < Mr < 15) in order to understand better 
their internal dynamics. Our goals are the extension of the 
Faber- Jackson relation (L oc cr") and the fundamental plane 
(FP, Djorgovski & Davis 1987; Dressier et al. 1987; Bender 
et al. 1992) relating luminosity, velocity dispersion, surface 
brightness, and scale length of galaxies, to the low luminos- 
ity dwarf population, providing important tools for studying 
formation of galaxies. For example, the precise value of the 
exponent of the Faber- Jackson relation, a, helps to constrain 
galaxy formation models as it is sensitive to the effect of 
gas loss from self-gravitating systems (Dekel & Silk 1986; 
Yoshii & Arimoto 1987). In the case where a dark matter 
halo governs the galaxy dynamics, a is predicted to be 5.26 
and M/L oc L^''''^ . If, instead, the galaxies contain only 
baryonic matter and have roughly constant M/L, a flatter 
relation, L oc u^ '^, is expected (Dekel & Silk 1986). 

Clusters of galaxies are ideal for studying scaling rela- 
tions as they provide large number of galaxies of different 
types, at a common distance, different environments within 
the cluster, and allow one to assess the interaction with the 
Intra-Cluster Medium (ICM). Since the Coma cluster is one 
of the nearest, rich and dense clusters, it has been a popular 
target for studies of scaling relations among galaxies includ- 
ing the Faber- Jackson relation (J0rgensen et al. 1996; Moore 
et al. (2002: MLKC02); Graham & Guzman 2003). 

The Faber- Jackson relation for faint ellipticals in Coma 
cluster has been studied by Matkovic & Guzman (2005: 
MG05) and Cody et al. (2009: Co09), using WIYN/HYDRA 
(multi-fiber) spectroscopy. Both MG05 and Co09 find Lr oc 
^2.o±o.3 jtqj. (j^ajf ellipticals even though these studies have 
only 10 galaxies in common whose a measurements differ 
by up to ~30% for the faintest galaxies. Thus, it is important 
to obtain high S/N data to measure velocity dispersions for 
dwarf galaxies at fainter magnitudes. 

In this paper we use the high spectral resolution of 
DEIMOS (Faber et al. 2003) on the Keck n telescope to 
measure the velocity dispersions for a sample of 41 faint 

elliptical galaxies in the core of the Coma cluster. Of these 
galaxies, 15 are common with the samples of MG05, Co09 
or both, and for the remaining 26, these are the first veloc- 
ity dispersion measurements. The current DEIMOS mea- 
surements extend the study to luminosities in the range 
— 16.5 < Mr < —15, about 1 magnitude fainter than the 
Umit reached by the WIYN/HYDRA data in previous stud- 
ies of Coma dwarf galaxies. 

This paper covers mainly the observations and the tech- 
nical part of the analysis. The emphasis is on the reliable es- 
timate of the uncertainties in measuring the internal velocity 
dispersion of galaxies. In section 2 we describe the observa- 
tions, the spectroscopic setup and the data reduction. Section 
3 is dedicated to the measurement of radial velocities and 
the velocity dispersions. The error analysis is described in 



section 4. The results of the analysis and the Faber- Jackson 
relation are presented and discussed in section 5. 

Throughout this paper, the distance modulus of the 
Coma cluster is considered to be 35.00 mag (Carter et al. 
2008). 



2 OBSERVATIONS 

2.1 Source Selection 

Galaxies were selected using the photometry from Adami 
et al. (2006). The (B-R) versus R colour-magnitude relation 
was constructed, and the red sequence identified and fitted 
by a linear relation. Galaxies were then selected from within 
a region ±3(7 from this relation, with a magnitude range 
-19 < Mr < -16, or 16 < R < 19. Images from the 
HST/ACS Coma Treasury Survey (Carter et al. 2008) were 
then used to reject a very small number of galaxies with ob- 
vious spiral structure. However not all candidates fell within 
the footprint of the completed survey (which was truncated 
owing to ACS failure) so some late-type galaxies may re- 
main in the sample. The masks were then designed at the 
university of Hawaii to accommodate, and as a compromise 
between the two programmes: the velocity dispersions pro- 
graimne described here, and a separate programme to mea- 
sure redshifts of a sample of low surface brightness galaxies. 
The requirement that the spectra from different slits should 
not overlap, and the need for sufficient slit length to enable 
sky subtraction, places considerable restriction on the num- 
ber of galaxies one can observe in each mask. 

2.2 Spectroscopic setup 

The observations were carried out on the night of March 
19"' 2007, using DEIMOS on the Keck II telescope. 
The outer dimension of the rectangular DEIMOS field is 
16.7x5.0 arcmin^ and all observed spectra are mapped onto 
2x4 mosaic of 2Kx4K CCDs (2 CCDs for each spec- 
trum). We used 1200G grating (1200 grooves/mm BK7, 
gold-coated, and blazed at 7760 A). The nominal central 
wavelength at the centre of the mask was 8700 A, although 
of course the central wavelength of a given spectrum de- 
pends upon the position of that slitlet in the mask. Although 
slitlets for the redshift survey observations were drilled at 
a variety of widths, all of those for the velocity dispersion 
project described here had widths of 0.7 arcsec. Recipro- 
cal dispersion is 0.33 A/pixel, giving a spectral resolution 
(FWHM) of 1.6 A, and a resolving power R~5000. This 
equates to an instrumental velocity width of ~25 km s~^, 
allowing us to measure velocity dispersions down to ~15 
km s~^ in spectra with S/N ^ 15 per pixel. The velocity 
dispersions of galaxies in the 17 ^ ^ 19 dE sample are 
expected to range between 10 km s~^ and 80 km s~^. The 
spectral range covered is approximately 7500 — 10, 000 A, 
although it depends upon the slit position in the mask. This 
wavelength range is dominated by the very deep calcium 
II triplet absorption lines (CaT) at 8498.02, 8542.09 and 
8662.14 A, which depend only weakly upon metallicity, and 
therefore provide an ideal range over which to measure the 
kinematic properties of low-luminosity galaxies. A GG495 
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Figure 1. Example of galaxy spectra in the sample with different values of signal-to-noise ratio (S/N). All S/N values are estimated per pixel. 
The Godwin, Metcalfe & Peach (1983; GMP) ID and R magnitude are given for each galaxy. In each panel, the arrow shows the location of 
CCD gap in the spectrum. 



filter was used to block out the contamination from the sec- 
ond order blue spectrum. 

Two masks were observed for each of two pointings, 
and the exposure time for each mask was in the range 6030 
to 6600 seconds. A summary of these observations is pre- 
sented in Table [T] In addition, a number of potential tem- 
plate stars, and flux calibration stars, were observed using 
a 0.7 arcsec width long slit during twilight at the beginning 
and end of the night. Wavelength calibration arcs and obser- 
vations of internal flatfielding lamps were also made at the 
beginning of the night. 

We observed eleven standard stellar spectra of which 
only two templates were reliable for measurements in this 
study. The other stellar templates were excluded because 
they either were not G/K-type stars or were affected by noise 
near the Ca triplet lines. 



pipeline. For our purpose here, all spectra are not flux cali- 
brated and all partially remaining sky emission lines, partic- 
ularly near the Ca triplet lines, were manually removed. In 
this manual process, for each galaxy spectrum, we checked 
the wavelength of the remaining bright emission lines with 
the catalogues of the sky emission lines and then we replaced 
the verified sky emission lines with the continuum fitted un- 
derneath the spectrum. Each spectrum was originally located 
on 2 CCDs. In order to join the blue and red parts of each 
ID spectrum, the average value of the nearest 50 pixels to 
the CCDs gap was calculated and the offset was used to 
construct the conjoined thorough ID spectrum. An exam- 
ple of final one-dimensional spectra from the DEIMOS in 
the range ~8300-9000 A is shown in Figure[I] In this paper 
all spectral S/N values are estimated per pixel (i.e. 0.33 A). 



2.3 Data reduction 

The spectra of dwarf galaxies and template stars were re- 
duced using the IDL spec2d pipeline developed by the 
DEEP2 Galaxy Redshift Survey team at the University 
of California-Berkeley for that project (see Davis et al. 
(2003) for more information). The pipeline performs the flat- 
fielding using internal quartz flats and the wavelength cali- 
bration using the ArRrNeXe arc lamps. The sky subtraction 
and the data extraction to ID spectra are also done within the 



3 THE MEASUREMENT OF RADIAL VELOCITY 
AND VELOCITY DISPERSION 

There are 49 galaxies in our DEIMOS sample with spec- 
troscopically confirmed Coma membership (Marzke et al. 
2011). The CaT absorption lines for eight galaxies were ei- 
ther coincident with gaps in the CCD mosaic, or too low 
signal-to-noise to make reliable measurements and therefore 
were taken out of the sample. We thus, used 41 spectra with 
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Figure 2. An Example of pPXF fit to a galaxy (GMP3424) from the 2007 DEIMOS run. Thin black is the observed spectrum, and the thick 
black (red) is the best pPXF fit using a single stellar template. The residuals over the fitted region is also shown (blue). The S/N for this spectrum 
is 29, and we find K = 5596 ± 7km s"!; cr = 27 ± 9km s~^. 



Table 1. List of the exposures on the galaxy masks. In column 5, 
we give the number of slits drilled at 0.7 arcsec width, specifically 
for the velocity dispersion programme, excluding those drilled at 
larger widths for redshift measurements of low surface brightness 
galaxies. In column 6, we give the total number of slits per each 
mask. 



Mask 


Mask Centre (J2000) 


Exposure 


No. of Slits 




RA 


Dec 


Time (s) 


dEs 


Total 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


Comal -1 


12:59:57.83 


27:59:30.6 


6600 


24 


105 


Coma 1-2 


12:59:57.83 


27:59:31.2 


6600 


15 


109 


Coma2-l 


13:00:04.04 


27:54:52.2 


6600 


22 


101 


Coma2-2 


13:00:04.04 


27:54:52.2 


6030 


14 


103 



identified CaT lines with signal-to-noise ratio (S/N) ranging 
from 5 to 87. Figure [T] shows four of the sample spectra. 

The line-of-sight radial velocities, and velocity disper- 
sions, a, were measured using the pPXlS software devel- 
oped by Cappellari & Emsellem (2004). The pPXF tech- 
nique works in pixel space and uses Gauss-Hermite series 
to extract the radial velocity and velocity dispersion simul- 
taneously by minimizing the which is defined as 

n — l ' ' 

where G(xn) and Gmod[xn) are the original and the 
modelled galaxy spectra respectively. AG(a;„) is the error 



^ penalized PiXel Fitting 



in observed galaxy spectrum and N is the number of good 
pixels used in the fitting process. 

First, both galaxy and stellar template spectra are re- 
binned in wavelength space (a; = ln[\)). A range of model 
spectra are obtained by convolving the stellar template with 
a broadening function, which takes into account the true 
galaxy radial velocity and a range of trial velocity disper- 
sions. The minimisation is then carried out over this range of 
broadened templates. Compared to the other methods, Cap- 
pellari & Emsellem showed that for low S/N, the use of max- 
imum penalized likelihood suppresses the noise effect in the 
solution. This makes the pPXF method robust for objects 
with low S/N such as our dwarf galaxies. 

pPXF enables us to use any desired wavelength range 
for measurements (see Figure [2] for example). In order to 
benefit from high S/N part of the spectrum, we chose the rest 
frame spectral range of 8450-8700 A which covers the CaT 
features. All galaxy spectra and stellar templates are contin- 
uum subtracted. Due to the noise, for a few faint galaxies 
{Mr > — 16 and S/N < 15), the measured velocity dis- 
persion is very sensitive to the selected wavelength range 
and therefore, high errors are expected in the final results. 

In order to provide realistic initial values of radial ve- 
locity to pPXF, the IRAF FXCOR package was used to ob- 
tain the relative velocity for each selected pair of the galaxy 
and template spectra. FXCOR is based on the Fourier cross- 
correlation technique developed by Tonry & Davis (1979). 
The location of the highest peak of the correlation func- 
tion represents the most probable blue/red-shift value of the 
galaxy. Although the FXCOR task can measure velocity dis- 
persions, it also requires simulations of multiple spectra in 
order to construct a proper relation between the width of the 
highest peak of the correlation function and the real velocity 



© 2010 RAS, MNRAS 000, 1 



Coma Dwarf Galaxies: I. Velocity Dispersion Measurements 5 



Table 2. The list of stellar templates, their spectral type and their 
stellar parameters (i.e. effective temperature, metallicity and surface 
gravity). T^y f and g are in terms of [° A'] and [cm/s^]. All stellar 
parameters are taken from Indo-US library. These templates were 
used to perform the initial measurements and simulations. 



Library 


lU [rlUJ 


or i irt, 




[re/rlj 


log(g) 


DEIMOS 


44007 


G5IV 


4850 


-1.71 


2.00 


(This Study) 


45282 


GOV 


5280 


-1.52 


3.12 




44007 


G5IV 


4850 


-1.71 


2.00 




45282 


GOV 


5280 


-1.52 


3.12 




117876 


G8III 


4582 


-0.50 


2.25 




161797 


G5IV 


5411 


0.16 


3.87 


Indo-US 


223-65 


G2IV 


5350 


-0.59 


3.50 




142198 


KOIII 


4700 


-0.08 


2.99 




48433 


KlIII 


4425 


-0.24 


1.35 




6497 


K2III 


4421 


0.02 


2.80 




78479 


K3III 


4509 


0.57 


2.54 



dispersion of the galaxy. Such simulations are not necessary 
when using pPXF. Other advantages of the pPXF package 
include the use of multiple stellar templates, where these are 
optimized and weighted, and shorter computing time. 

Of the eleven stellar templates observed by DEIMOS, 
only two spectra (i.e. HD44007 and HD45282) had high 
enough S/N near the strong CaT absorption lines to be used 
in our analysis. Therefore, we added more G/K-type spectra 
from the Indo-U.S. Coude Feed spectral library (Valdes et 
al, 2004). For practical reasons, especially computing time, 
we used 1 1 stellar templates (see Table|2} to perform single 
template analysis and simulations. To find the best optimal 
synthetic template, a set of 80 mixed templates were used to 
finalize the results. 

The first pPXF run shows a significant difference in 
the velocity dispersion values when different templates in 
Table |2] are used. The scatter in a values gives an esti- 
mate of the uncertainty due to the multistellarity of galax- 
ies. To reduce the effect of stellar population on the mea- 
sured sigma, a set of 50 mixed stellar templates (hereafter, 
Set50) from Indo-US library consisting of G/K-stars were 
simultaneously given to pPXF and the best combination was 
obtained by optimizing the values. For each galaxy spec- 
trum, the Set50 pPXF run significantly improves upon a sin- 
gle template, the values are improved by 10% (on aver- 
age), compared with 7% intrinsic statistical deviation of a 
single template pPXF run. This also points at the multistellar 
nature of galaxy spectra. The impact of the stellar template 
mismatch is discussed in ^4.2\ In 44% (20%) of the cases, 
the is improved by more than 10% (20%). 

To explore the sensitivity of the final results to the num- 
ber of templates, a set of 80 templates (hereafter Set80), con- 
sisting of Set50 and 30 more G/K-type templates, also cho- 
sen from Indo-US library, were used. Comparing the result- 
ing velocity dispersion for Set50 and Set80 indicates that 
for 36 of 41 galaxies, the discrepancy in the results is less 
than 2% and it does not exceed 10% for the rest. Although a 
larger number of templates improves the results, it dramat- 
ically increases the computation time. Moreover, the use of 
Set80 instead of Set50 does not improve the significantly 
(i.e. ~1% on average), and thus we did not continue adding 
more templates to SetSO. We find that all galaxies are mod- 



Table 3. The Hst of most frequent stellar templates used in the anal- 
ysis with Set50 and SetSO. These templates are used to reproduce 
more than two galaxy spectra. Columns 3, 4 and 5 show the effec- 
tive temperature, metallicity and surface gravity of each star respec- 
tively. Stellar parameters are taken from Indo-US library. The Hst is 
sorted in terms of the metallicity. Columns 6 and 7 indicate the uti- 
lization frequency of the templates of Set50 and SetSO. Set50 does 
not contain 7 of the templates available in SetSO. 



ID 


SPTYPE 


Te// 


[Fe/H] 


log(g) 


S50 


S80 


[HD] 




[K] 




[cm/s^] 


(%) 


(%) 


63791 


GOI 


4700 


-1.81 


1.70 


17 


10 


44007 


G5IV 


4850 


-1.71 


2.00 


5 


7 


1S7111 


G8IV 


4429 


-1.54 


1.15 


24 


22 


45282 


GOV 


5280 


-1.52 


3.12 


24 


29 


81192 


G8III 


4582 


-0.70 


2.75 


7 


2 


210752 


GOV 


5910 


-0.64 


4.25 


12 


20 


G.241-18 


G5V 


5511 


-0.61 


4.00 


- 


10 


G.165-11 


GOV 


5725 


-0.56 


4.00 


20 


17 


201099 


GOV 


5872 


-0.50 


4.06 


- 


7 


157214 


GOV 


5676 


-0.41 


4.33 




7 


99167 


K5III 


3930 


-0.38 


1.61 


10 


10 


82210 


G5III 


5250 


-0.34 


3.42 


7 


5 


837S7 


K5III 


4000 


-0.18 


1.60 


17 


22 


169414 


K2III 


4450 


-0.16 


2.67 


7 


7 


186486 


G8III 


4980 


-0.11 


3.08 


12 


15 


225212 




3700 


-0.03 


0.80 


10 


2 


150680 


GOIV 


5825 


0.00 


3.80 




7 


237903 




4070 


0.00 


4.70 


46 


20 


17925 


K2V 


5091 


0.10 


4.60 


10 


7 


161797 


G5IV 


5411 


0.16 


3.87 


10 


10 


63302 


K3I 


4500 


0.17 


0.20 




10 


182293 


K3IV 


4486 


0.25 


3.00 


20 


17 


G.196-9 


K5V 


4000 


0.28 


4.50 




24 


130705 


K4III 


4350 


0.51 


2.10 


22 


20 



elled with a combination of at least two templates. For 80% 
of the galaxies in our sample, more than two templates de- 
scribe the spectra and in 30% of the cases, a mixture of more 
than 4 templates were chosen. The results of the SetSO pPXF 
run are assigned to each galaxy and used for the rest of this 
paper (See Table [6j. 

When we used Set80 (SetSO), 38 (14) templates are not 
chosen to describe any of sample galaxies. 6 of 26 used tem- 
plates from Set50 are not chosen again by pPXF when using 
SetSO. Of 30 added templates to Set50, only 14 templates 
contribute to describe at least one galaxy spectrum in SetSO 
pPXF run. In table [3] a list of stellar templates which are 
used to reproduce more than two galaxy spectra is provided. 



4 ERRORS AND UNCERTAINTIES 

In this section we discuss the accuracy of our measurements 
and the sources of uncertainties. 

4.1 Statistical Uncertainties 

One important source of uncertainty is the statistical error 
which is highly correlated with the S/N. We adopt a "boot- 
strapping" technique to quantify this uncertainty. We sim- 
ulate an ensemble of galaxy spectra covering the observed 
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S: HD45282 (Deimos), M: HDl 17876 (In, US) 



S/N=10 
S/N=15 
S/N=30 
S/N=80 
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(a) 



S: HD45282 (In, US), M: HDl 17876 (In. US) 



3/N=10 
3/N=20 
S/N=30 
S/N=80 




40 50 60 70 80 90 100 110 120 130 
Input Velocity Dispersion:(7in (Km/s) 

(b) 



Figure 3. The RMS of measured velocity dispersion is plotted against the input velocity dispersion ((Ti„) used to run simulations. 300 simulated 
galaxies are generated for each S/N and (Ti„. In both plots, the templates are the same. In Plot (a), DEIMOS and Indo-US spectra used as 
simulation and measurement templates respectively. As one expects, lower S/N and greater (Ti„ result in bigger inaccuracies. In plot (b), both 
templates were selected from Indo-US library. Comparing the 2 plots, it is clear that when using different instruments for the simulation and 
measurement, the typical value of the uncertainties is greater. 



range of S/N, the broadening parameter, (cr), and the veloc- 
ity shift (Vr), using the stellar template spectra. To achieve a 
desired (ct), the stellar template spectrum is convolved with 
a Gaussian function. Then random noise is added to achieve 
the desired S/N. For each observed galaxy spectra, a set of 
300 realizations were generated, noting that, increasing the 
number of realizations does not alter the results. We then 
apply the same measurement technique. 

The measurements on each simulated set are carried out 
with the same initial parameters used in the measurement 
of original spectrum. Although the input spectral broaden- 
ing, (7i,i, for each set of simulated spectra is the same, the 
measured value ((Xm) can be different, due to the noise, and 
forms a Gaussian distribution around a mean value. The 
mean of the measured values is not necessarily the same as 
the input (7i„ value, due to a systematic shift. In general, 
for all templates in Table[2]and all probable S/N values, we 
obtained P{(Jm\o'in) (i-e. the probability of obtaining Om 
for simulated galaxy with the initial velocity dispersion of 
(Jin). To perform the simulations, for any given signal-to- 
noise ratio, (Tin is chosen to have a range of discrete val- 
ues 10-100 km . For other values of ain, we used the 
polynomial functions to obtain the mean and RMS of the 
measured velocity dispersion. This is explained in detail in 
section 3.3 of Co09. Figure [3] shows the relation between 
the input velocity dispersion and the RMS of the measured 
a for different S/N values. As expected, Om is more uncer- 
tain when the S/N decreases or velocity dispersion increases. 
Left panel of Figure [3] indicates that the uncertainty in cr is 
higher when the template for the velocity dispersion mea- 
surement is from a different instrument compared to that for 
the simulated galaxy spectra. 

In reality, we have only one spectrum of each galaxy 
based on which we estimate its velocity dispersion, am- 
Similarly, assume that of each set of simulated spectra with 
given S/N ratio and a in, we only have access to one spec- 
trum and its measured velocity dispersion, (7,„. Therefore, 
to derive the initial velocity dispersion (ain) and its un- 



certainty, we have to inverse the above statistics and find 
P{o-in\crm)- This is what needs to be done for any given 
original galaxy spectrum. Having all P((Jm\o'in,i) for each 
ain,i £ {1, 2, 3, 200 fcm s^^}, the probability of finding 
the true velocity dispersion, atrue for any measured velocity 
dispersion, Um, is given by (Co09) 

200 

P{o'tru<;,k\o'm) = -P(o"m 1 0"i„,fc ) / P((T,„ | Cri„,i ) , (2) 

i = l 

where (Jin,i ~ crtrue.i = i km s~^. Figure]?] shows 4 
examples of the probability density of true velocity disper- 
sion for any given measured value, a^- In order to illustrate 
the dependency of the results upon the selected templates, 
for panels "a" and "b" of Figure]?] we used the stellar tem- 
plate HD161797 (G-star) and for panels "c" and "d" we used 
HD78479 (K-star). As seen, the peak centre of the probabil- 
ity distributions are not centred upon the true velocity dis- 
persion. This systematic shift is mainly due to the template 
mismatch and other systematics caused by using different in- 
struments for simulation and measurement. As seen in Fig- 
ure]?] the location of the peak centres changes when we use 
different templates for measurements and for a given tem- 
plate it is independent of S/N. For any given velocity disper- 
sion and stellar template, the S/N value affects mainly the 
statistical uncertainty and therefore the width of the Gaus- 
sian peaks. 

To estimate the statistical inaccuracy of the final re- 
sults, la uncertainty is derived for each peak of the prob- 
ability density. The simulated galaxy sample is generated 
based on the same templates used for the measurements. 
We performed the above procedure for all individual stel- 
lar templates in Table ]2] which are initially used to measure 
a and V,- of the sample galaxies. Since all sample galax- 
ies are observed with DEIMOS, in another run of simula- 
tions we applied DEIMOS stellar templates (i.e HD44007 or 
HD45282) for simulation and Indo-US templates for mea- 
surement. This method makes our analysis more realistic 
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Figure 4. The distribution of true velocity dispersion for measured dispersion, am ■ Different instruments and templates are chosen to consider 
all sources of uncertainty. For each (Ji„ = 10, 20, 100 km s~^, 300 simulated spectra are generated using HD45282. In each diagram, the 
S/N of the simulated spectra is constant. The signal-to-noise ration in left and right panels are 10 and 30, respectively. To see the dependence 
of the uncertainties upon the stellar template used for measurements, measurements are carried out with HD161797 (G-star) and HD78479 
(K-star) in top and bottom panels, respectively and then the statistics are inverted to find P{crtrue\crm)- The width of the Gaussian peaks 
con'esponds to the statistical uncertainties (see il4.U and the systematic shift of the peaks along the axis of true velocity dispersion are due to 
the fact that we used different templates and different instruments for simulation and measurement (see i|4.2| for more detailed analysis). We 
see higher uncertainties for larger am values. To estimate the final uncertainty of the measured dispersions, la width of each corresponding 
peak and its systematic shift are added in quadrature. 



and enables us to study any additional sources of error due to 
the instrumental differences. When we apply individual stel- 
lar templates, the change in measured velocity dispersions 
for each galaxy is larger than the corresponding statistical 
error and we attribute this to the effect of the template mis- 
match which is investigated in ^4.2\ 



4.2 Template Mismatch and Other Source of 
Uncertainties 

The velocity dispersion measurement is highly sensitive to 
the width of the spectral lines in the stellar template based 
on which the measurements are carried out. The width of the 
lines is determined by chemical compositions, temperature 
and the age of the stars (Schulz et al. 2002; Lejeune at al. 
1997 and references therein). On the other hand, any galaxy 



is a composite of different stellar populations and thus lack 
of any prior knowledge about the exact dominant stellar type 
is an additional source of uncertainty. 

Velocity dispersion measurements of galaxies using the 
red part of the spectrum are reported to be less sensitive 
to systematic errors compared with measurements using the 
blue and visual spectral regions (Pritchet 1978). Our mea- 
surements are based on CaT lines which depend weakly on 
metallicity in metal rich stars (Mallik 1994). However, for 
metal poor stars, the effect of the metallicity on CaT lines 
could not be neglected (Starkenburg et al. 2010; Jorgensen et 
al. 1992; Zhou 1991; Smith & Drake 1990; Di'az et al. 1989) 
and therefore CaT lines are also employed as a metallicity 
indicator (Foster et al. 2010; Carrera et al. 2007; Michielsen 
et al. 2007). Moreover, evidence has been found for signif- 
icant anti-correlation between Ca II triplet indices and cen- 
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simulation by: HD 45282 (Delmos) Simulation by: HD 45282 (Indo-US) 
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Figure 5. The relative difference between the measured velocity dispersion and the given input value is plotted for different templates. In 
panel (a), DEIMOS spectrum HD45282 is used for simulation, and measurements were carried out with Indo-US and DEIMOS spectra. The 
template spectrum is broadened with the desired broadening factor (o"in) ™d no noise is added in this stage to see only the effect of template 
mismatch. Curve A shows the results of the same template used for simulation and measurement (i.e. DEIMOS HD45282) and as seen, the 
discrepancies are negligible. The measurements using another DEIMOS template (i.e. HD44007) leads us to curve B with larger discrepancies 
for low dispersions. Other curves show the measurement results using Indo-US templates. Panel (b) is the same as panel (a) except for the 
instrument used for simulation. As expected, there is no discrepancy when the same template is used for measurement (i.e. Indo-US HD45282). 
The black dash-dotted line in right panel shows the results of measurements with Deimos template (HD45282). As a result, the effect of 
mismatch and other uncertainties are less than 15% for > 40 km . Ignoring both templates HD44007 and HD45282 whose results 
behave differently regardless of the selected instrument, the average mismatch error is about 7%. We noted that the metallicity of these two 
templates are significantly less than the other templates (see Table|2). 



S: HD45282 (Delmos), M: HD151797 (In. US) 



S: HD 45282 (In. US), M: HD 161797 (Indo.US) 
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Figure 6. Same as the Figure |5]except that we used only one of the stellar spectra for velocity dispersion measurements and the S/N ratios 
of the simulated spectra vaiy. The relative differences of measured and input velocity dispersions are plotted for different S/N values. This 
shows that for S/N > 10 per pixel, mismatch and other uncertainties are not affected by signal level which means that the effect of S/N enters 
only through the statistical uncertainties. 300 simulated spectra were generated for each pair of S/N and (Ti,j . The simulations for left and right 
panels were carried out with DEIMOS and Indo-US spectra of HD45282. The same template (HD161797) was used for measurements in both 
panels. 



tral velocity dispersion of elliptical galaxies (Cenarro et al. 
2003). Therefore, template mismatch may have a significant 
impact on the final results of a measurements and, although, 
our measurements are based on 80 G/K-type stellar tem- 
plates, we still need to quantify the uncertainty due to the 
template mismatch. 

Another source of uncertainty appears when the galaxy 
and template spectra are observed with different instru- 
ments. This error is due to the difference in the sensitiv- 



ity and resolution of the instruments. The full width at half 
maximum (FWHM) spectral resolution of the Indo-US tem- 
plates is 1.0 A and that of DEIMOS is about 1.6 A for the 
wavelength range in this study. To take this into account, a 
proper Gaussian function is convolved with the higher reso- 
lution spectrum corresponding to the difference in the spec- 
tral resolution. In addition, the sampling rate of Indo-US and 
DEIMOS spectra are 0.4 A/pixel and 0.3 A/pixel, respec- 
tively. Although pPXF re-samples and re-bins both template 
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[Fe/H] 

Figure 7. The difference between the measured velocity dispersion, 
cTm , and the given input velocity dispersion, (7i„ , against the metal- 
licity, [Fe/H], of the templates used for measurements. We used 
the star spectrum HD45282, observed by DEIMOS, as the input 
spectrum. In each panel, this spectrum was broadened according to 
the a range of input velocity dispersions. For all simulated spectra, 
S/N=30. In each case, the error bar is calculated based on 300 sim- 
ulations. For measurements, we used several templates from Indo- 
US library with the same spectral type (i.e. G8III) and metallicities 
ranging between -1 and 0.4 (in solar units). In each panel, RMS is 
the root-mean-square of the difference between the input and the 
measured velocity dispersions, and CC shows the coefficient of the 
correlation between Cm — o"in and the metallicity of the templates 
used for measurements. As seen, the measured velocity dispersions 
are independent of the metallicity within ~ 2 km uncertainty 
for > 10 km s~^. 



and target galaxy spectra to account for this difference, this 
still may lead to systematic errors in the final results. 

Due to the limited number of stellar templates observed 
with DEIMOS, we are unable to disentangle the uncertainty 
due to instrumental differences from the uncertainty due to 
template mismatch. As seen in Figure ID the cumulative ef- 
fect of both template mismatch and difference in the instru- 
ments (Indo-US and DEIMOS characteristics) causes sys- 
tematics in the measured values. Therefore, the shifts seen 
between the ai„ and am is due to both effects. Similar dia- 
grams are generated for any selected Indo-US and DEIMOS 
stellar templates. All pairs of different templates and instru- 
ments are considered to study all sources of error and their 
contributions. 

We generated a comprehensive set of simulated data us- 
ing all stellar templates in Table [2] according to the recipe 
explained in section 4. 1 . In this method, measurements and 
simulations are carried out with different templates of the 
same instrument or with the similar templates from differ- 
ent instruments. The resulting velocity dispersions and given 
initial values are, then, compared after performing the cross 
measurements on the simulated data with various noise lev- 
els. Two examples are presented in Figure|5]which show the 
deviation of the measured values, am, from the true input 
values, ain, for a given velocity dispersion. All simulations 
are carried out without adding any noise in order to explore 
the template mismatch effect. We obtained higher systematic 
errors by choosing the templates from different instruments. 
All curves in the left panel of Figure[5]illustrate this fact ex- 
cept for two curves A and B for which DEIMOS templates 
were used for measurements. 

As expected, no significant systematic shift arises when 
the same template from the same instrument are used for 
simulation and also the measurement. Typically, the rela- 
tive systematic errors tend to be constant and to slightly 
increase, respectively, for higher and lower velocity disper- 
sions. To study the sensitivity of the systematics to the S/N, 
all aforementioned simulations are repeated with random 
noise added. Left and right panels in Figure [6] show the re- 
sults for the same curves in corresponding panels of the Fig- 
ure |5] including different noise levels. As seen, except for 
very low S/N (i.e. S/N ^ 10 per pixel), the systematic er- 
rors are not affected by value of S/N. We repeated the proce- 
dure several times for different wavelength ranges and any 
pairs of templates and instruments in Table [2] In all cases, 
for a > 20 km the effective mean value of all uncer- 
tainties does not exceed 20% of the true value of the velocity 
dispersion. 

Two examples are presented in Tables|4]and|5] In Table 
|5] we noticed higher measured a values for HD44007 and 
HD45282 compared to the measured values using the other 
templates, even though HD44007 and HD161797 (i.e. the 
template used for simulation) have the same spectral type. 
We attribute this to lower metallicity (i.e. [Fe/H]) of both 
HD44007 and HD45282 compared to the other templates 
(see|2l and consequent systematic errors. As seen in the sec- 
ond and third columns of Table |5] (low metallicity stars with 
[Fe/H]< —1.5), for atn > 20 km s~^, the average rela- 
tive discrepancy between the measured values, am, and in- 
put values, ain, is 20%. This value is reduced to 10% for the 
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Figure 8. Same as the Figure|2]except the templates used for simu- 
lations and measurements. To perform the simulation, we used sev- 
eral templates from Indo-US library with the same spectral type 
(i.e. G81I1) and metallicities ranging between -1 and 0.4. In each 
panel, the templates were broadened according to the desired in- 
put velocity dispersions. In this level of simulations, no noise was 
added to the simulated spectra. For measurements, we used the 
SetSO templates. As seen, when using SetSO, the measured veloc- 
ity dispersions are the same as the input velocity dispersions within 
~ 2 km uncertainty for (Ti„ > 10 km s~^. In each panel, 
RMS is the root-mean-square of the difference between the input 
and the measured velocity dispersions, and CC shows the coeffi- 
cient of the correlation between am — cri„ and the metallicity of 
the templates used for measurements. No significant correlation is 
seen between the measured and input velocity dispersions when we 
use SetSO. To avoid any confusion, the error bars are not shown in 
this diagram. Here, the typical uncertainty in measured velocity dis- 
persions is ~ 2 km s~^. 
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Figure 9. Radial velocities from MG05 (open circles) and MLKC02 
(filled diamonds) and Co09 (red stars) versus velocities in this study. 
The dotted line shows equality indicating that the velocities of 
MLKC02 and Co09 are systematically greater. 



other templates. Further evidence for systematic uncertain- 
ties in the presence of metallicity differences is seen when 
measuring the velocity dispersion of our sample galaxies 
with HD44007 and HD45282 in a single template run which 
results in 25% larger values, on average, compared to those 
measured using multiple templates. This discrepancy is re- 
duced to ~5%, on average, for single pPXF run with the 
other 7 templates of Table[2] 

In order to examine the effect of the galaxy and stel- 
lar template metallicity on measured velocity dispersions, 
we used 21 stellar templates with the same spectral type 
(i.e. G8III) from Indo-US library and different metallicity. 
As seen in Figure |7] we used the stellar spectra HD45282, 
observed by DEIMOS, to simulate galaxies with desired ve- 
locity dispersions, i.e. atn = 10, 30, 50, 80, 100 km s~^. 
In each case, 300 simulated spectra with S/N=30 were cre- 
ated. Then, the velocity dispersion of the simulated galax- 
ies were measured using the stellar templates with the same 
spectral type but different metallicity. Figure, |7] shows that 
differences between the measured and input velocity dis- 
persion are independent of the template metallicity within 
~ 2 km s^^ uncertainty in measured velocity dispersions. 
Due to the template mismatch effect, in each panel of the 
Figurel?] there is an offset between the given velocity disper- 
sion and the mean value of the measured velocity dispersions 
using the templates with different metallicity. This offset is 
due to the template mismatch effect and for a given spectral 
type, it seems to be independent of the template metallicity. 

In Figure [S] we used the mixed stellar templates, i.e. 
SetSO, for measuring the velocity dispersion of simulated 
galaxies with different metallicities. To perform the simu- 
lations, we broadened the stellar templates with different 
metallicities (i.e. the G8III-type stars used in Figure [7]l ac- 
cording to desired input velocity dispersions. At very low 
velocity dispersion (i.e. ain = 10 km s~^) the mean 
value of the measured velocity dispersion for all simulated 
spectra with different metallicities is 11.5 km For 
cr > 10 km s~^, the measurements result in the same 
value as the input values and are independent of metallic- 
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Figure 10. The difference in measured velocity dispersion in tliis study and common galaxies in tlie literature (MG05, MLKC02 and Co09) 
are plotted against our measured velocity dispersion and S/N value in left and right panel respectively. To derive Icr errors, the errors of our 
(T-measurement and those from literature are added in quadrature. The horizontal axis is logarithmic scaled. As seen, except 4 measured values, 
the rest are consistent with literature within the error bars. There is a good consistency for (T > 50 km for which (cr/crnt) is0.97it0.05 
and is 0.97 ±0.17 for o- < 50 kms'^. 



Table 4. This table shows the measured velocity dispersion using different stellar templates with different spectral types (see the first row). The 
first column represents the input velocity dispersion based on which the simulation is carried out. Under each template, the measured values are 
reported in terms of km s~^. DEIMOS star spectrum HD44007 is used to generate the simulated galaxy spectra and all measuring templates 
are picked up from Indo-US library. For each set, the average of the results is calculated and no noise is added to only estimate the effect of 
template mismatch which is the offset between the initial and measured values. In the last column, for each row the maximum discrepancy (i.e 
Amai = \crm — crin\/<^in (%)) is Calculated for (Tin > 10 (fcm in which (Ti,j is the initial velocity dispersion used for the simulation 

and Um is the measured value. For ai„ > 20 {km s~^), the average discrepancy is about 7% and for ai„ > 100 {km s^^) it is almost 3%. 
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ity. As a conclusion, based on the simulations, using the 
mixed stellar templates results in more reliable measure- 
ments and reduces the effect of template mismatch. To con- 
firm these results, we measured the velocity dispersion of 
the simulated spectra in all panels of Figure |7] using SetSO. 
For panels of Figure |7] from top to bottom, the results are 
(j„ = 10.5 ± 2.3, 29.7 ± 1.3, 49.8 ± 1.3, 79.5 ± 1.9 and 
99.3 ± 2.5 km s~^, respectively. As expected, there is no 
significant offset from the input velocity dispersion when we 
use the mixed stellar templates. 



provided by pPXF (e.g red curve in Fig.|2ll. This technique 
effectively gives the statistical errors of both measured ra- 
dial velocity and velocity dispersion. In the case of hav- 
ing few templates (e.g. when using Set50 or SetSO), this 
method is applied with the use of the equivalent synthetic 
template based on the best output weights. Applying the 
multiple stellar templates reduces the effect of the template 
mismatch particularly when dealing with low S/N spectra of 
faint galaxies which are expected to have smaller mass and 
internal velocity dispersion. 



4.3 Best-Fit Error Analysis 

In addition to the above, we repeated the same analysis for 
each galaxy to construct an ensemble of simulated spectra 
by adding the corresponding random noise to the best fits 



4.4 The Catalogue Uncertainties 

All the stellar templates in Table[2]were used to measure ve- 
locity dispersions for each galaxy in our sample. We used the 
RMS scatter around the mean values as the likely template 
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Table 5. Same as the Table |4] The template HD161797 is used for simulation and all templates are selected from Indo-US library. Looking at 
column 5, as one expects, no discrepancy between the initial and measured value is found. The two last columns (Aa & A;,) show the maximum 
discrepancies as in Table |4] Column Aa is calculated for all applied templates and since the templates HD44007 and HD45282 result in large 
discrepancies, the column A;, is created ignoring them. Considering all templates, for (7i„ > 20 {km s~^), the average discrepancy is about 
10% and for > 100 {km s ~ ^ ) it is almost 7%. The metallicity of HD44007 and HD45282 are lower than the other templates and therefore 
results in measuring larger am compared to the other templates (see Table[2). 
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mismatch uncertainty (see column 9 of the Table |6}. Since 
the exact combination of the stellar templates is unknown for 
each galaxy, the exact value of the template mismatch uncer- 
tainty depends on the selected templates for measurements. 
Even though Table [2] covers various metallicities and spec- 
tral types, choosing another set of templates may slightly 
alter this estimate, due to the multistellarity of galaxies. For 
66% of the galaxies in Table |6l the mismatch error is esti- 
mated to be less than 20% and only exceeds 25% for 22% 
of the galaxies. The average template mismatch eiTors are 
4%, 12% and 18% for galaxies with the measured ve- 
locity dispersions greater than 100 km s^^ , between 50 & 
100 km s^^ and between 20 & 50 km s^^ , respectively. 
These results which are purely obtained from observation 
are consistent with those estimated using simulations de- 
scribed in i]4.2l Figure |5] shows that the mismatch effect de- 
creases as the velocity dispersion increases. Simulations in- 
dicate that the average uncertainty due to the template mis- 
match is about 15% for cr ~ 50 km s^^ and ~5% for 
<j ^ 100 km , which is consistent with the reported 
values in column 9 of the Table |6] 

In ij3] it is mentioned that using Set50 improves the 
by ~10% on average, while adding 30 more templates re- 
duces the by only ~1% on average. Hence, a larger num- 
ber of templates was not considered due to its very limited 
effect. Since the reported a is based on the SetSO and the 
corresponding statistical uncertainty (see i ]4.3l l is not very 
different from that of the individual template analysis (see 
i]4.1b , we obtained the statistical uncertainties for each sam- 
ple galaxy by applying the method described in i]4. 11 and 
i]4.3l for SetSO (see column 10 of the Table [SJ. Assuming 
that the use of the SetSO eliminates the mismatch effect (see 
i|4.2b , one can take these values as the final uncertainty of the 
measured velocity dispersions. Being more conservative, in 
column 8 of Table [6] we derived the final error of the mea- 
sured a as the quadratic combination of both uncertainties 
(columns 10 & 1 1). These values are used for further analy- 
sis in this paper and paper II (Kourkchi et al. 201 lb). 

We performed the same analysis to obtain the uncer- 



tainty of the measured radial velocities. The effect of the 
template mismatch on measured radial velocities is less than 
1%, because this measurement is only sensitive to the exact 
location of the absorption lines of the applied stellar tem- 
plates which are slightly altered from one template to an- 
other. The final radial velocities and their corresponding un- 
certainties are obtained based on the SetSO run (columns 5 
& 6 of the Table©. 



5 RESULTS AND DISCUSSION 

We present the radial velocities, Vr, and velocity disper- 
sions, (7, of 41 dwarf elliptical galaxies in the Coma cluster 
in Table|6] Of these galaxies, 12 have S/N < 15 per pixel 
in their spectra. The reported values are measured based on 
SetSO, as explained in ij3] 

We compare our measured Vr and a values with MG05, 
MLKC02 and Co09 in Figures |9] and [TO] The mean dif- 
ference between our Vr values and those of MG05 is ~S 
km s~^. However when we include the Vr measurements 
from Co09, the mean difference is ~20 km s^^. Compar- 
ison of the velocity dispersion with MG05 (see Fig. llOt 
shows that for smaller values of velocity dispersion (i.e. 
a < 50 km s~^), our values are relatively smaller. Such 
discrepancy is likely due to the different methods of mea- 
surement, and/or to different templates used for measure- 
ment and/or to different wavelength regimes used for analy- 
sis. Since CaT lines are more prominent in dE galaxies than 
the absorption lines in the wavelength range ^^4200-5700 A, 
as used in studies of MG05 and Co09, our measurements of 
Vr and a are more robust. Our analysis also includes mea- 
surements with 80 mixed stellar templates and detailed anal- 
ysis of the uncertainties. 

In addition, of 41 studied galaxies in this paper, 12 
galaxies were observed twice using two different masks. In 
Table|7] for each of these galaxies we compare the measured 
velocity dispersions using both observed spectra. We also 
added both spectra of these galaxies to get better signal-to- 
noise ratios and then repeated the measurements. As seen. 
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Table 6. Radial velocity, Vr and velocity dispersion, cr of dwarf galaxies in our sample. In columns 10 and 11, we have a list of common 
galaxies in literature and their cr. In column 12, the data denoted by asterisks are from MLKC02 and the rests are from Co09. The likely en'ors 
in measured cr due to the effect of template mismatch are illustrated in column 9. The statistical uncertainties in a are listed in column 10. 
Column 8 shows the final uncertainty of the cr for all uncertainty sources and is obtained by adding the columns 9 and 10 in quadrature. Column 
13 is the S/N of the galaxy spectra which are calculated per pixel (0.33 A). Column 14 shows the number of stellar templates used to describe 
the galaxy spectra in Set80 run. 



GMP 


RA 


Dec 


^sdss 


C2 


Scz 


a 


A<7 


6cri 


d(T2 


cr (AfGOS) 


cr {Lit.} 


S/N 


#of 


ID 


(J2000) 


(J2000) 


mag 


km 


km 


km s^^ 


km 


km s^^ 


km 


km s^^ 


km s^^ 


1/pix 


Temp. 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


3534 


12:59:21.40 


27:58:24.80 


16.13 


6001 


7 


30 


7 


7 


3 


53 ± 3 


65 ± 6* 


45 


7 


3474 


12:59:26.40 


27:59:18.30 


18.69 


7666 


6 


22 


4 


3 


2 






37 


5 


3471 


12:59:26.60 


27:59:54.50 


15.35 


6624 


13 


79 


5 


5 


2 


86 ± 3 


85 ± 8* 


80 


7 


3438 


12:59:28.50 


28:01:09.40 


18.05 


5952 


31 


26 


8 


7 


3 


32 ± 11 




16 


3 


3424 


12:59:29.30 


27:56:32.00 


18.27 


5596 


7 


27 


9 


7 


5 






29 


6 


3406 


12:59:30.30 


28:01:15.10 


17.41 


7116 


7 


36 


3 


2 


2 






29 


3 


3376 


12:59:32.10 


27:55:15.80 


17.25 


7004 


22 


37 


3 


3 


2 




28 ± 4 


19 


5 


3340 


12:59:35.20 


27:56:05.00 


17.73 


4476 


12 


43 


8 


6 


5 






17 


2 


3336 


12:59:35.50 


27:54:21.60 


17.61 


6954 


6 


42 


4 


4 


2 


52 ± 5 




47 


5 


3325 


12:59:36.00 


27:54:22.00 


17.58 


8953 


13 


38 


7 


6 


2 






18 


2 


3312 


12:59:37.00 


28:01:07.00 


17.42 


7156 


27 


31 


4 


3 


2 


85 ± 3 




44 


5 


3308 


12:59:37.20 


27:58:19.90 


18.09 


7522 


8 


50 


4 


3 


2 






57 


4 


3296 


12:59:37.90 


27:54:26.40 


14.3 


7948 


21 


178 


7 


5 


4 


179 ± 4 


180 ± 4* 


87 


1 


3292 


12:59:38.00 


28:00:03.70 


16.25 


4951 


8 


40 


5 


3 


3 


63 ± 6 




48 


5 


3223 


12:59:42.40 


28:01:58.60 


18.41 


7797 


16 


22 


9 


8 


5 






11 


3 


3209 


12:59:44.20 


28:00:47.00 


18.58 


7111 


17 


30 


5 


4 


2 


66 ± 14 




30 


5 


3166 


12:59:46.90 


27:59:30.90 


17.39 


8292 


7 


38 


6 


6 


3 


58 ± 8 


35 ± 4 


26 


4 


3146 


12:59:48.60 


27:58:58.00 


17.74 


5358 


10 


38 


11 


9 


6 






11 


3 


3141 


12:59:49.10 


27:58:33.90 


19.86 


4957 


26 


59 


15 


11 


10 






5 


2 


3131 


12:59:50.20 


27:54:45.50 


17.29 


7183 


14 


19 


7 


6 


4 






22 


4 


3119 


12:59:51.50 


27:59:35.50 


19.74 


7060 


14 


37 


10 


9 


5 






8 


3 


3098 


12:59:53.90 


27:58:13.70 


17.73 


6740 


8 


33 


7 


5 


5 






28 


5 


3080 


12:59:55.70 


27:55:03.80 


18.11 


6599 


10 


9 


8 


6 


5 




16 ± 4 


19 


5 


3018 


13:00:01.00 


27:59:29.60 


18.17 


7464 


10 


38 


5 


5 


3 






15 


3 


2983 


13:00:04.00 


28:00:30.70 


19.03 


6349 


10 


27 


6 


5 


3 






18 


3 


2960 


13:00:05.40 


28:01:28.00 


15.18 


5866 


7 


60 


5 


5 


3 


68 ± 2 


64 ± 5* 


62 


5 


2931 


13:00:07.10 


27:55:51.50 


17.59 


7704 


17 


35 


6 


4 


4 






47 


5 


2877 


13:00:11.40 


27:54:36.40 


18.71 


7225 


12 


29 


9 


8 


5 






14 


5 


2839 


13:00:14.70 


28:02:26.90 


14.43 


5706 


15 


168 


6 


5 


3 


160 ±3 




57 


4 


2808 


13:00:17.00 


27:54:16.10 


19.64 


8944 


16 


69 


12 


9 


8 






13 


2 


2780 


13:00:18.70 


27:55:12.70 


19.06 


6554 


18 


63 


12 


9 


8 






13 


2 


2755 


13:00:20.20 


27:59:37.60 


18.87 


6812 


13 


27 


9 


8 


4 






14 


2 


2736 


13:00:21.70 


27:53:54.80 


16.41 


4857 


5 


35 


4 


3 


2 




35 ±3 


42 


4 


2718 


13:00:22.70 


27:57:55.00 


18.76 


6360 


12 


30 


8 


7 


4 






13 


4 


2676 


13:00:26.20 


28:00:32.00 


17.86 


5516 


30 


37 


8 


7 


4 






19 


4 


2655 


13:00:27.90 


27:59:16.50 


19.46 


9412 


14 


45 


5 


5 


2 






12 


3 


2654 


13:00:28.00 


27:57:21.60 


14.87 


7009 


21 


144 


7 


6 


3 


142 ±4 




66 


3 


2605 


13:00:33.30 


27:58:49.40 


18.28 


5126 


12 


36 


10 


8 


6 






13 


3 


2591 


13:00:34.40 


27:56:05.00 


17.17 


8653 


14 


53 


9 


8 


4 






20 


2 


2571 


13:00:36.60 


27:55:52.20 


18.68 


5944 


18 


17 


6 


5 


3 






20 


3 


2563 


13:00:37.30 


27:54:41.10 


19.5 


5895 


11 


25 


10 


9 


5 






9 


2 



considering the statistical uncertainties, the measurements 
are consistent. 



5.1 Faber Jackson Relation 

The relation between luminosity and velocity dispersion of 
bright elliptical galaxies was originally discovered by Faber 
& Jackson (1976) and is expressed as L oc cr". For brighter 
ellipticals a is ~4 while fainter galaxies exhibit a shallower 
slope. Davies et al. (1983) were the first to note that a 
changes from ~4 to ~2.4, for elliptical galaxies fainter than 
AIr ~ 21.7 mag. This result was also confirmed by Held 
et al. (1992) who found a — 2.5 for dE galaxies. These re- 
sults were finally extended to lower luminosities by MG05 
and Co09 who found a ~ 2.0. In this paper, we present data 
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for 26 dE galaxies in the Coma cluster down to AIr ~ — 15 
mag. Using the orthogonal distance regression, we derive 
the F-J relation to be L oc (j2-34±o.i9 entire sample, and 
L oc |ji-99±o i4 j^Qj. galaxies brighter than Mr ~ — 16 mag. 

We show the L — a relation for the galaxies in our sam- 
ple in Figures [TTI and [T2l where we also include the sam- 
ple of Co09 and MG05. We excluded the spiral galaxies 
from Co09 when deriving the F-J relation. L — a correlation 
coefficient is -83%, at a significance of 99.99% (i.e., only 
a 0.01% chance that the correlation is due to the random 
scatter). We perform linear fits via two methods, the least 
square bisector fit and the orthogonal distance regression. 
Both methods incorporate the errors in our measurements 
and those reported in the literature. In orthogonal distance 
regression, the orthogonal distance of data points from the 
line is minimized. The magnitude of all galaxies in Figure 
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Figure 11. The Faber- Jackson relation for three samples of dwarf galaxies in the Coma cluster. Filled black diamonds represent dwarf elliptical 
sample from DEIMOS/Keck observations. Red asterisks are dEs from Cody et al. (2009) and green open circles are data from Matkovic & 
Guzman (2005). Except one point, the DEIMOS sample provides the only measurements at Mji > —16. The dashed line is the least-square 
orthogonal distance fit for DEIMOS-I-MG05 data while the dash-dotted line is the best orthogonal fit for all three samples. Of our sample, six 
of faintest galaxies are labelled with their GMP IDs. Dotted line shows the least-square orthogonal fir excluding the six galaxies fainter than 
Mn = — 16 mag. See Table[8]for detailed results. 



Table 7. The measured velocity dispersion of the galaxies which 
were observed with two different masks. Columns 2 through 5 are 
the mask IDs also presented in Table[T] For these galaxies, we added 
both spectra of two masks to get the larger signal-to-noise ratios. In 
last column, the measured velocity dispersion of the added spectra 
is presented. The velocity dispersion are in terms of km s~^ and the 
errors are based on the statistical uncertainties. As seen, the mea- 
sured velocity dispersions for each galaxy are consistent within the 
error bars. 



GMP 


Comal-I 


Coma 1-2 


Coma2-I 


Coma2-2 


Add 


3424 






27±7 


26±7 


27±5 


3223 


16±9 


29±8 






22±5 


3209 


32±3 


29±5 






30±2 


3119 


40±7 


32±7 






37±5 


3080 






8±9 


13±9 


9±5 


2983 


23±6 


28±5 






27±3 


2877 






32±8 


30±7 


29±5 


2808 






75±11 


7I±9 


69±8 


2755 


22±6 


31±7 






27±4 


2718 






33±5 


29±6 


30±4 


2605 


41±9 


32±8 






36±6 


2571 






17±4 


14±4 


17±3 



nn and the left panel of Figure [12] are derived from SDSS 
DR7. To transform SDSS ugriz magnitudes into Johnson- 
Cousins R magnitude, we used the mean value of the follow- 
ing relations presented by Lupton 200^: 

ii = r - 0.1837 * (ff-r)- 0.0971, (3) 



^ For various transformations between SDSS ugriz magnitudes to 
Johnson-Cousins U BVRI system see Jordi et al. (2006), Ivezic et 
al. (2007) and 

\http://www. sdss. org/dr6/algonthms/sdss UBVRITmnsform. html\ 



7? = r - 0.2936 * (r - j) - 0.1439. (4) 

The results of the least square fits are available in Ta- 
ble [S] We run the same analysis for the same sample using 
the CFHT Megacam i-band data. We find that our FJ rela- 
tion slope is consistent with those of Co09 and MG05 within 
the errors, for both linear fitting methods. The least square 
bisector method yields Lr oc (j^'^^^'^"'^®, while we derive 
Lh oc (j^ ''^*'^'^'"* via the orthogonal regression method. 
Using the CFHT i-band data, the best fits of FJ relation 
are consistent with those derived from SDSS data (see Ta- 
ble [8} . The behaviour at the faint end of the FJ relation is 
uncertain as many of our galaxies in this magnitude range 
{Mr > —16) have S/N < 15 per pixel. However our re- 
sults imply that the slope of the FJ relation may change to- 
wards higher velocity dispersions at this faint end. Excluding 
galaxies fainter than Mr ~ -16, FJ is Lr oc ct^ '-'^^'' ". 

Next, we use photometry from the HST/ACS Coma 
Treasury survey (Carter et al. 2008) to derive the FJ re- 
lation. Out of 41 galaxies in our sample, 28 have avail- 
able HST/ACS magnitudes, together with 43 galaxies from 
MG05 and Co09 samples. In Figure[T2]and Table[8l we show 
the best linear fits of the FJ relation using HST/ACS F814W 
(I-band) and F475W (g-band). We find that the FJ relation 
shows no discrepancy on colour. 

Following all above analysis, we noted that for one 
of our sample galaxies (GMP 3308), the SDSS magnitudes 
are almost one magnitude fainter than what is derived from 
CFHT and HST images. Therefore, we modified the SDSS 
R magnitude of this galaxy by constructing a linear trend 
between the CFHT i-band and SDSS R magnitudes of all 
sample galaxies. 
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Table 8. Least square fit to find L — cr relation. ( 1 ) The sample galaxies used in regression. (2) Surveys-Filters for mag- 
nitude of the sample galaxies. (3) The order of regression. For example, in the first row residuals in loga are minimized. 
(4,5) The slope and intercept of M = A loga + B and their Icr uncertainties. (6) a in the Faber-Jackson relation 
(L oc 0-"). We have used the data with 5/Af > 15 per Aof MG05 and all dEs of Co09. 



Sample Galaxies 


Filter 


Regression 


Slope (A) 


Intercept (B) 


a 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


DEIMOS+MG05 


SDSS-R 


loga 


-6.59 ± 0.67 


-6.72 ±2.17 


2.64 ±0.27 


DEIMOS+MG()5 


SDSS-R 


Mr 


-4.55 ± 0.38 


-10.45 ± 0.76 


1.81 ± 0.15 


DEIMOS-HMG05 


SDSS-R 


Bisector 


-5.57 ±0.55 


-8.59 ±1.63 


2.23 ± 0.22 






Ovthogonal 


— O.DZ ± U.U / 




z.Oi it yj.z i 


DEIMOS*-hMG05 


SDSS-R 


Orthogonal 


-5.41 ± 0.49 


-9.02 ±2.02 


2.16 ± 0.20 


DEIMOS+MG05-hCo09 


SDSS-R 


Bisector 


-4.85 ± 0.40 


-10.12 ± 1.31 


1.94 ±0.16 


DEIMOS-hMG05-hCo09 


SDSS-R 


Orthogonal 


-5.85 ± 0.49 


-8.37 ± 1.79 


2.34 ±0.19 


DEIMOS * -hMG05-hCo09 


SDSS-R 


Orthogonal 


-4.98 ± 0.35 


-10.00 ± 1.58 


1.99 ± 0.14 


DEIMOS+MG05+Co09 


CFHTLS-i 


Bisector 


-4.73 ± 0.38 


-10.33 ± 1.27 


1.89 ± 0.15 


DEIMOS-hMG05-hCo09 


CFHTLS-i 


Orthogonal 


-5.61 ± 0.45 


-8.81 ± 1.73 


2.25 ± 0.18 


DEIMOS * -hMG05-hCo09 


CFHTLS-i 


Orthogonal 


-4.82 ± 0.33 


-10.29 ± 1.54 


1.93 ± 0.13 


[DEIMOS * -i-MG()5-hCo()9J 


HST/ACS-F814W(I) 


Orthogonal 


-5.55 ± 0.63 


-8.90 ±2.49 


2.22 ± 0.25 


[DEIMOS*-hMG05-hCo09] 


HST/ACS-F475W (g) 


Orthogonal 


-5.33 ± 0.61 


-8.02 ±2.35 


2.13 ± 0.24 


[DEIMOS * -^MG05-^Co()9] 


SDSS-R 


Orthogonal 


-5.37 ±0.64 


-9.21 ±2.64 


2.15 ± 0.26 



All galaxies in [DEIMOS-l-MG05-l-Co09] sample have HST/ACS data. DEIMOS* represents the DEIMOS sample 
excluding galaxies fainter than Mr = — 16 mag. 




.4CS Mf814w (I) [Mag] AC.S Mp-ivsw (s) [Mag] SDHK Mr [Mag] 

Figure 12. Same as the Figure [TT] The Faber-Jackson relation for three samples of dwarf galaxies which lie in the footprint of the HST/ACS 
Coma Treasury survey (Carter et al. 2008). Left panel shows the FJ relation for this sample using magnitudes based on SDSS data. For left 
and middle panels, magnitudes are derived from the Treasury HST/ACS images in the F814W (I-band) and F475W (g-band) filters (Paper II). 
Dashed line is the least-square bisector fit and dash-dotted line is obtained using the orthogonal distance regression. Dotted line shows the 
orthogonal fit for galaxies brighter than Mr = — 16 mag. We see no colour dependence of the Faber-Jackson within the error bars. 



5.2 Colour-cr relation 

The relation between dynamical mass of galaxies and their 
stellar population can be studied using the colour-cr dia- 
gram. [VIG05 showed a well-defined colour-cr relation for 
faint early type galaxies brighter than AIr = —17.5 mag. 
In panel (A) of Figure [T3] we extend the relation to fainter 
dEs (Mr < — 15 mag) which are located at the blue end of 
the diagram and are more scattered around the the fitted line 
for brighter galaxies (i.e. for Mr > — 17 mag the relation is 
B - R = (0.23 ± 0.03)Zoffcr + (1.02 ± 0.23)). The scat- 
ter of dEs around the colour-cr relation could be due to their 
different formation history, age and metallicity. To calculate 
the colour of the galaxies, B-R, we used the SDSS data and 



Lupton (2005) transformation functions to derive the John- 
son B and R magnitudes. Panel (B) of Figure [T3] shows the 
same relation using 771475 — ni8i4 colour of the galaxies in 
the HST/ACS field. 

The dependence of the internal velocity dispersion of 
dEs and their central activity is investigated in panel (C) 
of Figure [13] Central colour of the galaxies, for which we 
have the HST/ACS images, are calculated within the aper- 
ture with diameter of 0.4" (i.e. 185 pc at the location of the 
Coma cluster). We found a linear trend between the velocity 
dispersion of dEs and their central colours, except for very 
faint galaxies which have bluer core, (771475 — 7)78i4) < 1-0. 
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Figure 13. Colour-c relation for three samples of dwarf ellipticals in the Coma cluster. Filled black diamonds represent dEs from our 
DEIMOS/Keck sample. Open diamonds represent faintest dEs of our sample which are labelled with their GIVIP IDs. Red asterisks are dEs 
from Co09 and green open circles are from MG05. Dashed line shows the best Hnear orthogonal fit for brighter galaxies, Mjt > — 17 mag and 
dash-dotted lines represent \a uncertainty along the vertical axis. In all panels, 'a' and 'b' are slope and intercept of the fitted lines, respec- 
tively. Panel (A) shows the C-ir relation using the magnitudes from SDSS DRV. Panel (B) shows the C-cr relation for a sub-sample of dEs for 
which their HST/ACS magnitudes in F475W (g) and F814W (I) bands are available. Panel (C) is the same as panel (B) except for the colour of 
galaxies which are calculated within the aperture with diameter of 0.4" (i.e. 185 pc at the location of the Coma cluster). 



Table 9. The results of a measurements for 8 blue galaxies based 
on different template libraries. 



GMP 




a {km s ^ ) 




ID 


Set80 


FGK-lib 


F-hb 


3146 


38 


37 


35 


3141 


59 


60 


61 


3119 


37 


35 


26 


2877 


29 


27 


26 


2808 


69 


39 


38 


2780 


63 


63 


31 


2655 


45 


45 


42 


2563 


25 


25 


28 



5.3 Remarks on faint galaxies 

Of six faintest galaxies from Deimos sample, GMP2808, 
GMP2655, GMP3119, GMP2780 and GMP3141 are bluer 
than all other sample galaxies (see Figure 113b and have 
larger velocity dispersion with respect to the L-ct relation 
of brighter galaxies (see Figure[8l(. 

In order to test the effect of stellar population on mea- 
sured velocity dispersions, we have repeated the a measure- 
ment using two sets of stellar templates. First set consists 
of 10 F-stars with different spectral types and metallicities 
(hereafter F-lib). Second set includes both F-lib and the G/K- 
stars of Table [3] (hereafter FGK-lib). The results of a mea- 
surements based on SetSO, FGK-lib and F-lib are compared 
in Table[9] The analysis based on FGK-lib shows no changes 
in results except for GMP2808 which is the bluest galaxy in 
our sample. The measured a of GMP31 19 and GMP2808 in 
F-lib run is smaller than those of Set80 and FGK-lib runs, 
however the resulting is higher for F-lib run. The maxi- 
mum change in measured a for all other galaxies in our sam- 



ple in FGK-lib run is less than 2 km . As a results, the 
measured a values in Table |6] are independent of the stel- 
lar population except for the bluest galaxy, for which blue 
stellar templates should be used for a measurement. 

Recent gas-rich mergers in blue dEs may be responsi- 
ble for their deviation from the FJ relation. Supernova driven 
winds can expel the luminous matter of these galaxies and 
shift them to the faint end of the FJ diagram. Moreover, the 
resulting starburst due to the inward gas transportation dur- 
ing the merger, produces the central excess light compared 
to the extrapolated Sersic function of the galaxy outer region 
(Kormendy et al. 2009). The correlation between the central 
extra light and deviation of faint dEs from the scaling rela- 
tions of brighter ellipticals is studied in paper II. 



6 SUMIMARY 

We have observed a sample of ^^^50 dwarf elliptical galaxies 
in the core of Coma cluster of which for 41 dwarfs we mea- 
sured their internal velocity dispersion. For 26 galaxies, we 
presented the velocity dispersion for the first time. Our cur- 
rent study extends the relation one magnitude fainter than 
previous studies in Coma. We performed a comprehensive 
analysis to find the source of the uncertainties when mea- 
suring the velocity dispersion. We found that the main un- 
certainty may arise due to the mismatch of the stellar tem- 
plates. This kind of uncertainty is less than 20% of the mea- 
sured velocity dispersion for a 50 km s^^. The tem- 
plate mismatch uncertainty decreases for higher velocity dis- 
persions and is reduced to ~ 5% for a ^ 100 km s~^. 
We find greater velocity dispersions when using templates 
with lower stellar metallicities. The effect of template mis- 



© 2010 RAS, MNRAS 000, 1 



Coma Dwarf Galaxies: I. Velocity Dispersion Measurements 17 



match is reduced with the use of multiple stellar templates 
and the optimization of the template weights to reproduce 
each sample galaxy spectrum. In addition, using the multi- 
ple stellar templates removes the dependency of the mea- 
surements on galaxy and/or templates metallicity. Casting 
our data points and those from the literature, we get L-cr re- 
lation as L oc (ji-ssio 14 noticed that fainter dwarfs 
show a departure from the FJ relation of the brighter ellip- 
ticals which indicates that they have higher velocity disper- 
sion than what is predicted by L-cr linear trend. In paper II, 
we will present the combined study of the photometric and 
spectroscopic properties of this sample, focusing on the fun- 
damental and the photometric plane of the galaxies. 
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